Stable vortex-antivortex molecules in mesoscopic superconducting triangles 
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A thermodynamically stable vortex-antivortex pattern has been revealed in mesoscopic type I su- 
perconducting triangles, contrary to type II superconductors where similar patterns are unstable. 
The stable vortex-antivortex "molecule" appears due to the interplay between two factors: a re- 
pulsive vortex-antivortex interaction in type I superconductors and the vortex confinement in the 
triangle. 
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Symmetrically-confined vortex matter in superconduc- 
tors, super-fluids and Bose-Einstein condensates offers 
unique possibilities to study the interplay between the 
Coo symmetry of the magnetic field and the discrete 
symmetry of the boundary conditions. More specifi- 
cally, superconductivity in mesoscopic equilateral trian- 
gles, squares etc. in the presence of a magnetic field 
nucleates by conserving the imposed symmetry (C3, C4) 
of the boundary conditions 0] and the applied vortic- 
ity. As a result, in an equilateral triangle, for ex- 
ample, in an applied magnetic field H generating two 
flux quanta, 2$ , superconductivity appears as the C3- 
symmetric combination 3$o — (further on denoted 
as "3 — 1") of three vortices and one antivortex in the 
center. These symmetry-induced antivortices can be im- 
portant not only for superconductors but also for sym- 
metrically confined superfluids and Bose-Einstein con- 
densates. Since the order parameter patterns reported 
in Refs. Q have been obtained in the framework of the 
linearized Ginzburg-Landau (GL) theory, this approach 
is valid only close to the nucleation line T C (H). Can then 
these novel symmetry- induced vortex-antivortex patterns 
survive deep in the superconducting state? Several at- 
tempts have been already made to answer this crucial 
question. In the limit of an extreme type II supercon- 
ductor (k> 1), it has been shown that for a thin-film 
square, a configuration of one antivortex in the center 
and four vortices on the diagonals of the square is un- 
stable away from the phase boundary According to 
the analysis based on the coupled nonlinear GL equa- 
tions, the vortex-antivortex pairs are unstable and no 
antivortices appear spontaneously at the T, H points far 
away from the T C (H) line ||. Possible scenarios of pen- 
etration of a vortex into a mesoscopic superconducting 
triangle with increasing magnetic field have been stud- 
ied in Ref. Q]. Two different states were considered: a 
single vortex state and a state in the form of a symmet- 
ric combination of three vortices and an antivortex with 
vorticity L av = —2 ("3 — 2" combination). The calcula- 
tions B have shown that while a single vortex enters the 



triangle through a midpoint of one side, the "3 — 2" com- 
bination turns out to be energetically favorable when the 
vortices are close to the center of the triangle. Equilib- 
rium is achieved when a single vortex is in the center of 
the triangle. When approaching the phase boundary, the 
free energy of a single- vortex state tends to coincide with 
the free energy of the "3 — 2" combination Q , thus con- 
firming conclusions that formation of antivortices is 
possible in the close vicinity to the phase boundary. 

The previous inferences ^,|) on vortex-antivortex 
states in mesoscopic structures seem to give us no hope 
to find stable vortex-antivortex configurations deeper in 
the superconducting state, considering them just as the 
features appearing in materials with n 1 at the phase 
boundary together with superconductivity. Here we pro- 
pose the new solution demonstrating the stability of the 
vortex-antivortex patterns. This solution is based on the 
simple conjecture made by one of the authors (VVM, Q): 
the main source of the vortex-antivortex pattern instabil- 
ity, namely vortex-antivortex attraction, can be removed 
by taking - instead of type II - type I superconductors, 
where vortex-antivortex interaction becomes repulsive. 
Indeed, when passing through the dual point k = 1/V%, 
the vortex-vortex interaction changes the sign |^-|| and 
becomes attractive at k < 1/V2. At the same time, the 
vortex-antivortex interaction becomes repulsive. There- 
fore, one can expect that presence of antivortices, to- 
gether with confinement of vortices and antivortices due 
to a potential barrier at the boundaries, can stabilize 
novel vortex-antivortex patterns in a mesoscopic sample 
of type I superconductor. Optimizing the geometry and 
the sizes of mesoscopic samples, one can therefore ful- 
fil the conditions necessary for the existence of stable 
vortex-antivortex configurations. For instance, presence 
of sharp corners is known (^-jl^l to lead to a strongly in- 
homogeneous distribution of the superconducting order 
parameter in a mesoscopic sample. Enhanced supercon- 
ducting condensate density at the corners prevents vor- 
tices from leaving the sample. Altogether, a triangular 
type I superconducting sample seems to be an appro- 
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priate candidate to search for a stable vortex-antivortex 
configuration. 
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FIG. 1. Schematics of unstable (type II) and stable 
(type I) vortex-antivortex molecules in a mesoscopic super- 
conducting triangle. 

Fig. 1 further illustrates this idea. In a type I su- 
perconductor with £ > A, vortex cores overlap hrst when 
vortices approach each other. This triggers vortex-vortex 
attraction (and vortex-antivortex repulsion) . In contrast 
with that, in a type II superconductor with £ < A, 
local fields created by vortices overlap first when vor- 
tices approach each other, thus inducing vortex-vortex 
repulsion (and vortex-antivortex attraction). As a re- 
sult, a vortex-antivortex combination, consisting of three 
vortices (L^ v = 3) and one antivortex with vorticity 
L av = —1 ("3 — 1" combination, or 3w + law molecule) 
is stable or unstable, respectively, in a mesoscopic type I 
or type II superconducting triangle. 

To verify these intuitive considerations we investigate 
a mesoscopic equilateral triangle of a type I supercon- 
ductor in an applied magnetic field. The side of the tri- 
angle is chosen (a — 1/im) to be larger than both char- 
acteristic lengths of the GL theory, £ and A. In order 
to provide a sample to be a type I superconductor, we 
should consider in fact a triangular prism with a height 
h £ (and also h A). This allows us to avoid an 
increase of n resulting in a known effect when a meso- 
scopic sample made of material, which is type I super- 
conductor in bulk, becomes effectively type II supercon- 
ductor [0Jl5-|l3]. For comparison, a mesoscopic triangle 
of type II superconductor will be also considered. In our 
calculations, we use the GL parameters of Pb (type I): 
£pb = 82 nm, Apt = 39 nm, np^ = 0.48, and of Nb (type 
H) : ^Nb = 39 nm, A^ = 50 nm, « Nb = 1.28 fll5| . 

In the description of the superconducting properties of 



mesoscopic triangles we rely upon the GL equations for 
the order parameter ip and the vector potential A of the 
magnetic field H = rot A |l4| In the dimensionless 
form, when keeping the temperature dependence explic- 
itly, the GL equations are 



(— iV — A) ip — tp 



k 2 AA = - (V>*VV> - ipVip* 



= 0, 
A | VI 2 - 



(1) 
(2) 



Here the GL parameter k = A(T)/£(T), and £(0) serves 
as the unit length. The imposed boundary condition is 



n • (-iV - A) V| 



boundary 



= 0. 



(3) 



Topological characteristics of solutions of the GL equa- 
tions are determined by (anti)vortex core lines. One rev- 
olution along any closed path around such a line changes 
the phase of the order parameter by 27rL, where L is the 
winding number (vorticity) of a vortex or antivortex. 

Close to the sample, the local magnetic field is dis- 
torted as compared to the applied magnetic field H. 
However, far away from the sample the distortion is neg- 
ligible, and the symmetric gauge of the vector potential 
can be applied: A = [H x r] /2. The boundary condition, 
which equates the vector potential of the local field with 
that of the applied field, is justified at the boundaries 
of the simulation region, which is taken large enough to 
provide that all the changes of the magnetic field occur 
inside this region (see, for example, Refs. [[|jl7],[l|] ) . Fur- 
ther on, we choose the z-axis coinciding with the direc- 
tion of the applied magnetic field, which is normal to the 
triangular base. Because fe> A, a non-zero z-component 
of the magnetic field is appreciable in the vicinity of the 
bases, where the field skirts the prism partially penetrat- 
ing it. Near the central cross-section in the rry-plane, 
the magnetic field and the order parameter are uniform 
in the z-direction, and the problem becomes effectively 
two-dimensional. The third dimension is taken into con- 
sideration by imposing that, inside the simulation region, 
the total magnetic flux through any cross-section, per- 
pendicular to the z-axis, is the same 

• Thc GL 

equations (111), (||) with the above boundary conditions are 
solved numerically, using the finite-difference method, on 
a square mesh with the density of 200 nodes per side of 
the triangle. 

The calculations are performed for temperature 
T/T c — 0.96. First, we find the magnetic field regions, 
in which states with antivortices are expected to take 
place in a type I superconducting triangle. The free en- 
ergy calculations show that at low fields, from Hq = 
to 0.14iJ c (0), where H c (0) is the thermodynamical crit- 
ical field |14]Jl5|l at zero temperature, no vortices pen- 
etrate the triangle. The order parameter distribution 
changes continuously from homogeneous at Hq = to 
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a strongly inhomogeneous one, characterized by concen- 
tration of the superconducting phase in the corners. At 
Ho ~ 0.14ff c (0) one vortex enters the triangle. For vor- 
ticity L = 1, the equilibrium is achieved when the vortex 
is in the center of the triangle. For magnetic fields within 
the region from Ho — 0.25 to 0.38i/ c (0), a giant vortex 
state with the total vorticity L — 2 is energetically prefer- 
able. Further on, we focus on H = 0.32iJ c (0), which 
corresponds to the states with the total vorticity L = 2. 
This state can be represented by two possible configura- 
tions: (i) two vortices in the form of a multivortex or a 
giant vortex state; (ii) a symmetric combination of three 
vortices on the triangle bisectors and one antivortex in 
the center ("3 — 1" combination, or 3v + lav molecule). 
(Symmetric combinations with a larger number of vor- 
tices and antivortices such as "6 — 4", "9 — 7" etc. possess 
a higher free energy than the "3 — 1" combination and 
are not considered.) 
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distance d, from the center of the triangle, ^(0) 

FIG. 2. The free energy F s — F n [measured in Hi (0)/47r] as 
a function of the distance d v from the center of a mesoscopic 
type I superconducting triangle for a symmetric combination 
of three vortices and one antivortex (the "3 — 1" combination, 
or 3v + lav molecule), at T/T c = 0.96, H = 0.32# c (0), 
k = 0.48. 

According to our calculations, it is the "3 — 1" combi- 
nation that minimizes the free energy in case of a type I 
superconductor. In Fig. 2, the free energy for this combi- 
nation is shown as a function of the distance d v counted 
from the center of the triangle along the bisectors to vor- 
tices. There are three minima of the free energy as a 
function of d v . The first minimum, which is at —1.29^(0) 
from the center of the triangle, corresponds to a config- 
uration when vortices are situated between the center of 
the triangle and the midpoints of the sides of the triangle. 
This is a saddle point for the free energy as a function 
of the coordinates (a;, y) in the plane of the triangle, and 
the state is unstable. The second minimum is reached 
when all the vortices are in the center of the triangle, 



and the vortex-antivortex combination degenerates to a 
giant vortex Li v — 2. This local minimum represents a 
metastable state. 
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FIG. 3. The distribution of the squared modulus of the 
order parameter | if)(x,y) | 2 for the states with the total 
vorticity L = 2 in mesoscopic superconducting triangles at 
T/T c = 0.96, H = 0.32# c (0): 3u + lav molecule in a type 
I superconducting triangle with k = 0.48 (a); a giant vortex 
state in a type II superconducting triangle ft = 1.28 (b); a 
stable multivortex state in a type II superconducting triangle 
(c). 

The absolute minimum is reached when three vortices 
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are situated between the center and the apexes of the 
triangle at 1.83£(0) from the antivortex in the center 
(Fig. 2). This vortex- antivortex molecule is thermody- 
namically stable. Its stability can be understood in the 
following way. The distribution of the squared modu- 
lus of the order parameter | ip(x,y) | 2 , which relates to 
the above stable vortex-antivortex molecule, is shown in 
Fig. 3a. Four zeros of | ip(x,y) | 2 correspond to three 
vortices and one antivortex. The distribution of the mag- 
netic field H(x,y) consistent with the above pattern of 
| ip(x,y) | 2 will be presented elswhere pOj. The function 
| ip( x iy) | 2 reaches its maximum value in the corners. 
These "islands" of the superconducting phase in the cor- 
ners prohibit vortices, which are repelled by an antivor- 
tex in the center, from leaving the triangle through the 
corners. Thus, vortices, being confined in a type I super- 
conducting triangle and interacting with an antivortex, 
form a stable vortex-antivortex molecule. 

It is worth noting, that for our type I sample a strongly 
enhanced nucleation field appears due to the confinement 
of the superconducting condensate in the mesoscopic tri- 
angle. In fact, this provides the "soft" scenario for the 
nucleation of the order parameter, like in bulk type II 
superconductors. 
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FIG. 4. The free energy F s - F n [measured in H^(0)/4tv] 
for the state with the total vorticity L = 2 as a function 
of d v for a mesoscopic type II superconducting triangle at 
T/T c = 0.96, H = 0.32tf c (0), k = 1.28: 3u + lav molecule 
(a); a two- vortex state with vortices situated at two different 
bisectors of the triangle (b) . 

The stable vortex-antivortex patterns are qualitatively 
different in case of a type II superconducting triangle. 
In Fig. 4, the free energy of the "3 — 1" combination is 
shown (curve "a" ) as a function of the distance d v from 
the center of the triangle to vortices. The lowest mini- 
mum is reached when all the vortices are in the center 
of the triangle. This means that a giant vortex with 
vorticity L2 V — 2 is energetically more favorable in a 



type II superconducting triangle. (The local minimum at 
d v = 11.57£(0) represents an unstable state.) The corre- 
sponding distribution of the squared modulus of the order 
parameter | ip(x,y) | 2 is plotted in Fig. 3b. Although a 
giant vortex state with Li v — 2 has a lower free energy 
than the "3 — 1" combination, the equilibrium is reached 
for a multivortex state when two vortices are at two dif- 
ferent bisectors (cf. Ref. ||) of the triangle (see Fig. 4, 
curve "b"). The order parameter pattern corresponding 
to this stable two-vortex state is shown in Fig. 3c. 

In conclusion, we have found deep in the superconduct- 
ing state a thermodynamically stable vortex-antivortex 
configuration for a mesoscopic type I superconducting 
triangle, although until now it has been thought that 
vortex-antivortex patterns are unstable and they can 
manifest themselves only in the close vicinity to the phase 
boundary. Vortex-antivortex arrays become unstable in 
a type II superconducting triangle, in accordance with 
previous reports. The stability of the vortex-antivortex 
molecules in type I superconducting triangles is due to 
the change of the sign in the vortex-vortex and vortex- 
antivortex interaction forces when passing through the 
dual point k = l/\/2, combined with the condensate con- 
finement by the boundaries of the mesoscopic triangle. 
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